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Abstract 

We investigate microwave background anisotropies in the presence of pri- 
mordial magnetic fields. We show that a homogeneous field with fixed direction 
ON \ can amplify vector perturbations. We calculate the correlations of 5T/T ex- 

plicitly and show that a large scale coherent field induces correlations between 
&i-i,m an d Q>t+l,m- We discuss constraints on the amplitude and spectrum of a 

ON \ primordial magnetic field imposed by observations of CMB anisotropies. 

OO 
O 
!>• 

: 1 Introduction 

OO 

Since the detection of CMB anisotropies by the COBE satellite 0, it has become clear 
that anisotropies in the cosmic microwave background (CMB) provide a powerful tool 
for distinguishing models of cosmological structure formation. Furthermore, they may 
help to determine cosmological parameters which influence their spectrum in a well 
defined, non-trivial way 0. It is thus important to calculate the CMB anisotropies 
for a given model simply and reliably. 

The origin of observed galactic magnetic fields of the order of /xGauss is still 
unknown. For some time it has been believed that tiny seed fields can be amplified 
by a non-linear galactic dynamo mechanism. The effectiveness of this process has 
recently been strongly questioned, however ||. If magnetic fields are not substantially 
amplified by non-linear effects, but have just contracted with the cosmic plasma during 
galaxy formation, primordial fields of the order 10~ 9 Gauss on Mega-parsec scales are 
required to induce the observed galactic fields. 

It is interesting to note that a field strength of 10~ 8 Gauss provides an energy 
density of Qb — B 2 /(87cp c ) ~ 10~ 5 fi 7 , where f2 7 is the density parameter in photons. 
We naively expect a field of this amplitude to induce perturbations in the CMB on 
the order of 10~ 5 , which is just the level of the observed anisotropies. This leads us 
to investigate the extent to which the isotropy of the CMB may constrain primordial 
magnetic fields. It is clear from our order of magnitude estimate that we shall never 
be able to constrain tiny seed fields on the order of 10~ 13 Gauss with this approach, 
but primordial fields of 10~ 9 to 10~ 8 Gauss may well have left observable traces in the 
microwave background. 
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This is the question we investigate quantitatively in this paper. Some work on 
the influence of primordial magnetic fields on CMB anisotropies has already been 
published, particularly the cases of fast magneto-sonic waves and the gravitational 
effects of a constant magnetic field || [7], Here we study Alfven waves. We leave 
aside the problem of generation of primordial magnetic fields. This issue is addressed 
e.g. in |,§. 

The possible influence of magnetic fields on large scale structure formation has 
been investigated in || and references therein. 

The paper is organized as follows. In Section 2 we discuss the effect of a homo- 
geneous magnetic field background on the cosmic plasma. Both scalar (potential, or 
fast and slow magneto-sonic waves) and vorticity (Alfven) waves can be induced. We 
study the latter. In Section 3 we consider the influence of Alfven waves on CMB 
anisotropies. The influence of the magneto-sonic waves can be interpreted as a slight 
change in the speed of sound, and has been investigated in In Section 4 we present 
our conclusions. Some of the more technical computations as well as a discussion of 
Silk damping of vector perturbations are left to two appendices. 

Notation: For simplicity we concentrate on the case Qq = 1. The choice of Q 
is of little importance for our perturbation variables (which have to be calculated 
for at early times, when curvature effects are not significant), but does influence the 
resulting Ces due to projection effects. Throughout, we use conformal time which we 
denote by t. The unperturbed metric is thus given by ds 2 = a 2 (t)(—dt 2 + 5ijdx l dx^). 
Greek indices run from to 3, Latin ones from 1 to 3. We denote spatial (3d) vectors 
with bold face symbols. 

2 Cosmological vector perturbations and Alfven 

waves 

Vector perturbations of the geometry are of the form 

<*-">= U *j?jf J • « 

where B and H are divergence-free, 3d vector fields supposed to vanish at infinity. 
Studying the behaviour of these variables under infinitesimal coordinate transforma- 
tions (called gauge transformations in the context of linearized gravity), one finds 
that the combination 

cr = H — B (2) 

is gauge invariant. Geometrically, cr determines the vector contribution to the per- 



turbation of the extrinsic curvature [H], |Tl|. 



To investigate perturbations of the energy-momentum tensor, we consider a baryon, 
radiation and cold dark matter (CDM) universe for which anisotropic stresses are 
negligible. The only vector perturbation in the energy-momentum tensor is thus a 
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perturbation of the energy flux, u, the time-like eigenvector of T u ' . We parameterize 
a vector perturbation of u with a divergence free vector field v, such that 

u = -v . (3) 

a 

Analyzing the gauge transformation properties of v, one finds two simple gauge- 



invariant combinations [10 



V = v - H and fi = v - B . (4) 
They are simply related by 

V = n - a . (5) 

The perturbations of the Einstein equations, together with energy-momentum con- 
servation, yield |T 



--Act = 3(-) 2 Cl, (6) 



2 \a 
* + 2 ('£)*■ = 0, (7) 



n + (i-3c;)-n = o. 



a 



The two Einstein equations (|5], [7]) and the momentum conservation equation (Rp are 
not independent. Eq. (§) follows from Eqs. (||) and (0). 

This system does clearly not describe waves. From Eq. (|7|) it follows that cr decays 
like 1/a 2 . Furthermore, Eq. implies ft oc (kt) 2 cr. In the radiation dominated era, 
where a oc £, this yields O = \(kti n ) 2 <T in , where U n is some initial time at which 
fluctuations were created, e.g. the end of inflation. The fact that O stays constant 
during the radiation dominated era also follows from Eq. (§). On cosmologically 
interesting scales, k <C l/£m, we have therefore f2 <C crj n and cr <Tj n . In contrast, 
scalar and tensor perturbations remain constant on super-horizon scales. For this 
reason, vector perturbations induced at a very early epoch (e.g. inflation) which 
have evolved freely can be entirely neglected in comparison to their scalar and tensor 
counterparts. 

This situation is altered in the presence of a primordial magnetic field, which 
induces vorticity waves after the inflationary era. Let us consider a homogeneous 
magnetic field B before the time of decoupling of matter and radiation. Such a field 
could have originated, for example, at the electroweak phase transition ||. When 
the photon-baryon fluid is taken to be a perfectly conducting plasma, an external 
magnetic field induces two distinct modes of oscillation. Magneto-sonic waves, scalar 
perturbations which propagate at speeds c± slightly above or very much below the 
ordinary speed of sound in the plasma. These induce density oscillations just like 
ordinary acoustic waves. The slight change in the speed of sound can change slightly 
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the position and shape of the acoustic peaks in the CMB spectrum 0. Here we 
discuss the vectorial Alfven waves. 

We assume a plasma with infinite conductivity and use the frozen-in condition: 

E + v x B = , 

where v is the plasma velocity field. Our plasma is non-relativistic (v <C 1). The field 
lines of a homogeneous background magnetic field in a Friedmann universe are just 
conformally diluted, such that Bq oc 1/a 2 . Until recombination, the photon-baryon 
plasma is dominated by photons, p r ~ p 7 oc a~ 4 (p r denotes the combined baryon 
and photon energy density) and the ratio Bl/(p r +p r ) is time-independent. We study 
purely vortical waves which induce a fluid vorticity field f2(k) normal to k. We note 
that charged particles are tightly coupled to the radiation fluid and obey the equation 
of state p r = p r /3. 

It is convenient to rescale physical quantities like the fields and the current density 
as follows: 

E — > Ea 2 , B -»• Ba 2 and J — > Ja 3 . 

We now introduce first-order vector perturbations in the magnetic field (Bi) and in 
the fluid velocity (fi); 

B = B + Bi, V ■ Bi = and (9) 

v = 0, V-ft = 0. (10) 

To obtain the equations of motion for Q, and B 1; we first consider Maxwell's equations. 
Since the fluid velocity is small, we may neglect the displacement current in Ampere's 
law, which then yields 

J = — VxBl (11) 

47T 

We replace Ei with Bo, using the frozen-in condition. The induction law then gives 

l-Bi = V x (v x B ) . (12) 
at 



Inserting relation (11) for the current, the equation of motion (T 1 ^ = F^j^) for vector 
perturbations in the plasma becomes 

|-v = - 1 B x (V x Bi). (13) 
at 47r(p r +p r ) 

(We have neglected viscosity, which is a good approximation on scales much larger 
than the Silk damping scale ||12|| .) Taking the time derivative of this equation, we 
obtain with the help of Eqs. @, ( |10|) and (|T2|) for a fixed Fourier mode k 

0= (B °' k)2 n and (14) 
47r(p r +p r ) 
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• zBn ■ k „ . . 

ft = — — -Bi . 15 

These equations^ describe waves propagating at the velocity tu( e • k), where 

B 2 



v 2 = — , d a ~4x 10- 4 (So/10 _9 Gauss) 

4ir(p r +p r ) 



is the Alfven velocity and e is the unit vector in the direction of the magnetic field. 
Typically the Alfven velocity will be very much smaller than the speed of acoustic 
oscillations in the radiation-dominated plasma (c 2 = 1/3 ^> v\). 

Due to the observed isotropy of the CMB, we have to constrain the magnetic field 
contribution to the total energy density. For example, in the radiation dominated era 
it must be a fraction of less than about 10~ 5 [^], leading to va ~ 10~ 3 . Eq. fll4] ) is 
homogeneous in Q and so does not determine the amplitude of the induced vorticity. 
The general solution contains two modes, cos(t> Aktfi) and sin(t> A^tfi) (where /z = e-k). 
If the cosine mode is present, it dominates on the relevant scales k < 1/ (tu^dec)- Then 
we can approximate cos (vAkt dec fi) ~ 1 and the sine is negligibly small. But this mode 
then describes the usual vector perturbations without a magnetic field. We assume 
it to be absent. We want to consider initial conditions, then, with 

ft(k,t = 0) = 0. 

Only the sine mode is present and we have 

fi(k, t) = fl sm(vAktfi) ~ Q, VAktfi . (16) 

The initial amplitude of flo is connected with the amplitude of Bi by means of 
Eq. (Jl5|), yielding 

IfioHK/AOIBxl . (17) 
This allows a vorticity amplitude of up to the order of the Alven velocity (see also 



3 CMB anisotropics from Alfven waves 

Vector perturbations induce anisotropics in the CMB via a Doppler effect and an 



integrated Sachs- Wolfe term |TT 



(^V" C) = _V.„Ko 



/ Jt dec 



where the subscripts dec and denote the decoupling epoch (z rfec ~ 1100) and today 
respectively. Since the geometric perturbation cr is decaying, the integrated term is 

: Our derivation is valid either in a gauge invariant framework as outlined in |l3}| or in a gauge with 
vanishing shift vector. In other gauges metric coefficients will enter and complicate the equations. 
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dominated by its lower boundary and just cancels cr in V = f2 — a. Neglecting a 
possible dipole contribution from vector perturbations today, we obtain 

ST 

— (n, k) ~ n • n(k, tdec) = n • n v A kt dec (e ■ k). (19) 

We assume that the vector perturbations Qq are created by some isotropic random 
process, and so have a power spectrum of the form 

(O Qi (k)O 0i (k)) = (<% - • (20) 

For simplicity, we further assume that the spectrum A(k) = (l/2)|f2 | 2 (^) is a simple 
power law over the range of scales relevant here, 

k n 

A{k) = A —^, k < ko, (21) 
k 

for some dimensionless constant Aq and cutoff wavenumber ko- With this we can 
calculate the CMB anisotropy spectrum. 
The C/s are defined by 



'ST ST , 



= ^£(2^+1)QP|(m) • (22) 



A homogeneous magnetic field induces a preferred direction e and the correlation 
function (|22"D is no longer a function of n • n' = fi alone but depends also on the angles 
between n and B as well as n' and B . Statistical isotropy is broken. Setting 

ST 

-jr( n ) = E fl fa y fa( n ). (23) 

in the isotropic situation, the C/s of Eq. fl22|) are just 

C e = (a em a* e J , (24) 

where () denotes a theoretical expectation value over an ensemble of statically identi- 
cal universes. We find that the presence of the preferred direction Bo not only leads 
to an m-dependence of the correlators, but also induces correlations between the mul- 
tipole amplitudes a^ + i i?Ti and a^_i jm . Correlations in the temperature fluctuations at 
different points on the sky are no longer simply functions of their relative angular 
separation, but also of their orientations with respect to the external field. Detailed 
computations of the correlators for the Doppler contribution from Alfven waves are 
presented in Appendix A. We obtain 

Ci(m) = (a em a* em ) 



A v 2 a -p (Mo 



\t J v u u/ r(-n/2)'' 
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( 




(25) 



D e {m) 




n + l)/2) 2 



(£ + m + 1)(£ - m + 1)(£ + m){£ - m) 
(21- 1)(27+ l) 2 (2^ + 3) 



) 



1/2 



r(£ + n/2 + 3/2) 
r(^-n/2 + 1/2) ' 



(26) 



This result is valid in the range —7 < n < — 1. For n < —7 the quadrupole diverges 
at small k, and for n > — 1 the result is dominated by the upper cutoff ko, 



For n = —5 we obtain a scale-invariant Harrison-Zeldovich spectrum, Cg ~ £ 2 . 

To obtain some insight into the effect of the cross terms Dg, we picture the corre- 
lation function 



for various orientations of the magnetic field with respect to the fixed direction n . 
These are shown for the case n = —5 and with n = z in figures 1 to 3. Notice, 
however, that these figures do not represent temperature maps but are plots of the 
correlation function. For a given realization stochastic noise has to be added. The 
explicit expression for /(n) is given in Appendix A. 

With no a priori knowledge of the field direction, it could be inferred by performing 
CMB measurements with various n and comparing the obtained /(n) with the plots 
below. Of course this procedure suffers from problems with cosmic variance, as once 
we fix a direction no in the sky we have only a single realization with which to 
determine /. The expectation value in expression (|2"8|), then, strictly refers to a 
(hypothetical) average over an ensemble of universes. 

A probably simpler observational test of the existence of a constant magnetic field 
is the presence of temperature correlations for unequal I. To simplify, we introduce 
the mean values 




n > -1 . 



(27) 




(28) 




(29) 
(30) 



where the bar denotes average over different values of m, and we find 




for n < -1 (31) 
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Figure 1: An Aitoff projection of the function /(n) for a homogeneous magnetic field 
pointing in the 9 = 7r/4, <fi = it/ 2 direction and the reference vector no pointing in 
the z-direction {6 = 0) (see equation (|28|) ). 




Figure 3: The function /(n) for B pointing in the 8 = direction (i.e., parallel to 
no). The gray scale scheme has enhanced the variation in /. 



Ct * A (^) 2 (k Q t )- 2 v 2 A -±—£ 2 , farn>-l (32) 
V to / n + 1 

D t /Ut ~ 3/2. (33) 

The existence of significant correlations between the a^_i )m and a^+i, m is a clear in- 
dication of the presence of a preferred direction in the universe. Due to its spin-1 
nature, a long-range vector field induces transitions £ — > £ ± 1 and thus leads to the 
correlators Dg. 

There are no published limits on these cross correlation terms. Since the full a^ m 's 
are needed to obtain them, full sky coverage and high resolution, as will be provided 
by the MAP and PLANCK satellites, are most important to obtain such limits. The 
galaxy cut in the 4-year COBE data leads to an influence of Cg by Ce±2 which is 



on the order of 10%. for 2 < £ < 30 §L4 \. To be specific, let us assume that this 



be a limit on the off diagonal correlations Dg. Then, first of all, the observed CMB 
anisotropics are not due to Alfven waves, since 0.1 ~ Dg/Cg is substantially smaller 
than the figure in Eq. d33|) . To obtain a limit on the magnetic field amplitude and 
the spectral index, we now require 

£ 2 Dg < 0.1£ 2 Cg ~ 10~ n for 2 < £ ~ 100 . (34) 
We now argue as follows. From Eq. ( |T7|) , and the fact that Bi ~ B , we have 

\Q Q \ 2 k 3 £ v 2 A . (35) 
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This inequality must hold on all scales inside the horizon at decoupling, k > l/td ec - 
With Eq. (|2l| ) we therefore obtain 



2A (k/k ) n+3 £ v \ \/t dec < k < k Q , (36) 

which implies 

2A)(M dec )- (n+3) ~ v\ for n < -3, (37) 

2A ~ v\ for n > -3 . (38) 

Here we have identified k with the maximal frequency (cutoff) of the magnetic field, 
which has to be introduced in the case n > —3 for ft not to diverge at small scales. 
A definite upper limit on ko is the scale beyond which the magnetic field is damped 
away, due to the finite value of the conductivity. The physical damping scale is given 



by [15 



(k D /ay = Attct/t , (39) 

where r denotes the cosmic time (not comoving) and a is the plasma conductivity. 
The conductivity of a non-relativistic electron-proton plasma is easily shown to be 
a ~ 4T, and it has been shown recently that this result still holds approximately in 
the very early universe [1TB . 



Using Tdec ~ 0.3 eV ~ 0.6 x 10 4 cm 1 and Tdec ~ 10 5 years ~ 10 23 cm, we obtain 
the comoving damping scale at decoupling 



fcoftfec) ~ k D (t dec ) ~ (z^y^WTcTdJr^ ~ 3 x lO-^cm" 1 , (40) 

and 



(Mo)(W) = k (t dec )T /adec ~ T ^lQ7rT dec /T dec ~ 0.4 x 10 14 . (41) 
Inserting the limiting values of Eqs. ( ^71 P^ ) for the A in Eq. (^JJ), Eq. ( P4| ) yields 

S^^r 5 ^ 2 ^!;^ 1 ^ ^ 5 < 10- 11 for n<-3, (42) 

3^(2.5 x l0-^ +3) 3r( _ n ) 2)2 r+5 < 10 ' n fOT -3<n<-l, (43) 

3^^(2.5 x 10~ 14 ) 2 -^£ 4 < 10- 11 for -1 < n, (44) 

Using va ~ 4 x 10 _4 (I?o/10 -9 Gauss), this can be translated into a limit for B Q which 
depends on the spectral index n and the harmonic i. In Fig. 4 we plot the best limit 
on B as a function of the spectral index n. To optimize the limit we choose I = 2 for 
n < —5 and I = 100 for n > —5. For n > —3 the limit becomes very quickly entirely 
irrelevant due to the huge factor 10 14(n+3 ). This reflects the fact that for n > —3, the 
magnetic field fluctuations grow towards small scales, and B% ~ Bq is leads to a limit 
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-7 -6 -5 -4 -3 

n 

Figure 4: The upper limit on the magnetic field amplitude Bq due to CMB anisotropies 
caused by Alfven waves, shown as a function of the magnetic field spectral index n. 
We assume Di < 0.1CV Allowed values of the field must lie in the dashed region. 
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at the tiny scale k = ko; where as the CMB anisotropies, are caused by the smaller 
fluctuations at large scales, k ~ £/to- At n < — 1 the induced D/s start to feel the 
upper cutoff and thus do not decrease any further. 

The presence of a homogeneous magnetic field also induces anisotropic stresses 
in the metric. This gravitational effect has been estimated elsewhere J7|. Compared 
with the COBE DMR experiment, it leads to a similar limit for B . 

4 Conclusions 

We have studied Alfven waves in the primordial electron-proton plasma that are 
sourced by a homogeneous magnetic field. In addition, we allow for an isotropic dis- 
tribution of random magnetic fields on smaller scales. The induced vorticity in the 
baryon fluid leads, via the Doppler effect, to vector- type CMB anisotropies on all an- 
gular scales larger than the vectorial Silk damping scale (-damp ~ 500 (see Appendix B). 
The vector nature of the magnetic field induces off-diagonal correlations, 

D e (m) = (ai_i, m a| +1>m ) ~ C e (m) . (45) 

Assuming that observations constrain these terms to be less than about 10% of the 
observed C/s, we derive a limit for the amplitude of the magnetic field as a function 
of its spectral index. For n < — 7, the quadrupole anisotropy diverges if no lower 
cutoff is imposed on the spectrum, and so such spectra are very strongly constrained. 
For n > —3, the constraint is proportional to (&o£o) (n+3 ' ) ' /4 , where to is the comoving 
scale today and k is the upper cutoff of the spectrum. We have set l/k equal to 
the magneto-hydrodynamical damping scale which is inversely proportional to the 
conductivity and thus extremely small, leading to (k to) ~ 10 14 . Therefore, the limits 
obtained for n > —3 are extremely weak and actually uninteresting. This is due to 
the fact that the quantity B 2 (k)k 3 decreases on large scales for n > —3. For spectral 
indices in the range —7 < n < —3 the limit on B is on the order of (2 — 7) x 10 _9 Gauss. 

An important remark is also that causally induced magnetic fields lead to a spec- 
tral index n = 2 and so are not constrained at allQ. Examples here are magnetic fields 
generated by the decay of a Y field during the electroweak phase transition. 

At first sight it may seem somewhat artificial to have split the magnetic field 
into a homogeneous component and an isotropic spectrum of random magnetic fields. 
However, this is the correct procedure for realistic observations. This is seen as 
follows. If we calculate the Cgs for a given model, we determine expectation values 
over an ensemble of universes. If we make a measurement, however, we have just one 
observable universe at our disposition. This problem is generally referred to as 'cosmic 
variance'. On scales much smaller than the horizon, cosmic variance is irrelevant if 

2 If we assume magnetic fields to be generated by a causal procedure, i.e. not during an inflationary 
epoch, in a Friedmann universe, then the real space correlation function has to vanish at super- 
horizon distances (say \x\ > 2t). Its Fourier transform (BiBj){k) oc k n (5ij — kikj) is therefore 
analytic in k, which requires n to be an integer with n > 2. 
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we make some kind of ergodicity hypothesis, assuming that spatial and ensemble 
averaging are equivalent. In the case here cosmic variance is especially important, 
however. If we include the field which is coherent on the horizon scale in our random 
distribution, then in performing the ensemble average we integrate over the directions 
of Bo and the off-diagonal correlators vanish. The entire effect disappears. In one 
given universe, however, this field has one fixed direction and our effect is observable. 
It is therefore of fundamental importance here not to take an ensemble average over 
the large scale coherent field. The cosmic variance problem is relevant, whenever 
perturbations with non- vanishing power on horizon scales are present. 

An observational limit on the off diagonal correlators D/s from MAP or PLANCK 
would represent a model independent limit on the importance of large scale coher- 
ent vector fields (which enter quadratically in the energy momentum tensor) for the 
anisotropies in the cosmic microwave background. Its importance thus goes beyond 
the magnetic field case discussed in the present work. 
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A Calculation of Ci from vector perturbations 

In the usual way we decompose the temperature fluctuations of the microwave back- 
ground into spherical harmonics: 

Y"(n) = J2a irn Y em (n). 

Ira 

The two point function (a em a}, m ,) is then 

If r r /ST* ST \ 
(a em 4 m >) = J^jd ik j dttnj dSl* ^— (n, k) — (n', k) \ ^(n)y^(n') . (46) 

We consider the contribution to the temperature anisotropy only from the vorticity 
in the baryon fluid. From Eq. (|T9"D, we obtain 

^(n, k, At) = e 4k nA *n ■ ft(k, t dec ) , (47) 
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where At = to — td ec — h, the time elapsed since last scattering. Using the form ( p0| ) 
for the power spectrum of the vortical velocity fluctuations, we have 

^*(n,k)^(n',k)^ = ^< n '-^\n-n'-^')(v A kpt dec ) 2 A(k) 

= /(n,n',k), (4* 

where fj, — n • k, fj,' — n' • k, /3 — e • k and e is the unit vector in the direction of the 
homogeneous magnetic field B . So 

<<W*Jm'> = ^tfS d " k S dn " J dn "' n '' k ) ^m(n)^ m .(n') . (49) 
To evaluate these integrals, we use the identities 



4vr£ £ ^(*)y r ;(k)y r >) , (so) 

r=0 q=—r 



4tt +1 

n • n' = P 1 (n-n') = — £ ^(^(n') , (51) 

6 P =-i 

where j r is the spherical Bessel function of order r. Using the orthonormality of the 
spherical harmonics and the recursion relation 

(2i + l)f e = e Je -i-(£+l)je + i , (52) 



we find that in evaluating (|49"D only the terms with (£, m) = (£', m') and {£, m) 
(£' ± 2, m') survive, where 

/2£ 4 + 4£ 3 - £ 2 - 3£ + 6m 2 - 2£m 2 - 2£ 2 m 2 \ 
K a tm a tm) ~ y (2£-l)(2£+l) 2 (2£ + 3) J X 



and 



- [dkk 2 (v A kt dec ) 2 A(k)( Je+1 + j^) 2 (53) 

7T J 



(%+lm a f-l,m) ~~ ( a £-lm a £+l,m) ~ 

-(P-lVP-u* f (£ + m+l)(£-m + l)(£ + m)(£-m) \ 1/2 
1 ^ + ^ (2£-l)3(2£+l) 2 (2^ + 3) 3 J 

- [dkk 4 {v A t dec ) 2 A{k) ( Je + j £ _ 2 ) (j, + j e+2 ) . (54) 

7T J 
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The Bessel functions take kto as their arguments. With A(k) = Ao(k/ko) n k 3 , we 
obtain, for —7 < n < — 1, 

= C e (m) 

2 n+1 A v 2 A /t dec \ 2 T(-n - !)£(£ + n/2 + 3/2) 
(Mo)" +3 U / T(-n/2) 2 T(£ - n/2 + 1/2) X 



(2l 4 + 4l 3 - £ 2 - 31 + 6m 2 - 2£m 2 - 2l 2 m 2 ) 
(2£- l)(2£ + 3) ' 

( a t+l,m a t-l,m) = ( a t-l,m a t+l,m) = D t( m ) 



(55) 



2 n+2 A v\ ft dec \ 2 r(-n-l) T(£ + n/2 + 3/2) 
\n + l\(k t ) n+3 \t ) r(-(n + l)/2) 2 T(£ - n/2 + 1/2) X 

+ (2£-l)(2£ + l) 2 (2£ + 3) J • (56) 

For n > — 1 the integral is dominated by the upper cutoff k and we find 

Ce (m) = v A^o - C-^-Yx (57) 

V ; 2n(n + l)(k t ) 2 \t J K J 



{2£ A + 4£ 3 -£ 2 -3£ + 6m 2 - 2£m 2 - 2£ 2 m 2 ) 



(58) 



n / \ _ ^A^O (tdec\ 

UiKm) ~ 27r(n + l)(A;oto) 2 U / 



(2£-l)(2£ + 3) 

- I /+ \ 2 

x (59) 



^ (f + m + l)(£-m+l)(f + m)(f-m) V /2 
( ^ 1)(£ + 2) ( (2£-l)(2£+l) 2 (2£ + 3) J • (60) 

In this case, the result is nearly independent of the spectral index n and, due to the 
factor (kot )~ 2 , it is so small that it fails to lead to relevant constraints for B . 
The temperature correlation function is finally 

/ 5T ST \ 
/(n) = (— (n ) — (n)\ = J2 (a em a* e , m ,)Y em {ii )Y* m ,(n) 

Y,( a e+i,m a *e-i,m) (^+i,m( n o)^/-i, m (n) + V 1>m (n )F; +lim (n)) . 

im 
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B Collisional damping for vector perturbations 



Denoting the fractional perturbation in the radiation brightness by A4, Ai = 4(AT/T), 



the Boltzmann equation for vector perturbations gives [11 



M + n ■ VM = -AnWa^j +aa T n e [-M + 4n ■ il] . (61) 

Here n is the photon direction, a? denotes the Thomson cross section and f2 is 
the baryon vorticity. We have neglected the anisotropy of non-relativistic Compton 
scattering. 

To the baryon equation of motion (g) we have to add the photon drag force, 



a 3pb 

with 



-M- ft 

4 



(62) 



47T 

We shall also use the fact that for vector perturbations, the perturbation of the photon 
brightness vanishes, 

r Mdn = . 



3 

M = — / nMdn . 



Due to the loss of free electrons during recombination, the mean (conformal) collision 
time t c = l/(a<jTn e ) increases from a microscopically small scale before recombination 
to a super-horizon scale after recombination. After recombination the collision term 
can be neglected and we recover Eqs. (|3|) and (|IBfl. We first consider the very tight 
coupling regime, t c A, where A denotes the typical scale of fluctuations. In this 
limit the term inside the square brackets of Eqs. (BT) and ( [62] ) can be set to zero and 
we obtain M = 4f2 ( the baryon and photon fluids are adiabatically coupled). 

Next, we derive a dispersion relation for the damping of fluctuations due to the 
finite size of t c . We proceed in the same way as Peebles |17| for scalar perturbations. 



We consider scales with wavelength k 1 t and thus neglect the time dependence 
of the coefficients in Eqs. (|5"T|) and To study the damping we also neglect 



gravitational effects, which act on much slower timescales. With the ansatz 

M = A{n)exp(i{\sx-ut)), (63) 
CI = £(n)exp(z(kx-cjt)), (B ■ k) = (64) 



we obtain 



-iuA + i(kn)A = -[-A + An-B] (65) 
= ^|M-4B], R S %±. (66) 
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In the limit kt C} ujt c — > 0, we again obtain adiabatic coupling. The general relation 
between A and B is 

4n- B 

A = 

1 + i(k ■ n — uj)t c 

and so 



M = 3B 



(kt c 



-((kt c ) 2 + (1 - iut c ) 2 ) In L _ + + 2zfa c (l - iut c 



(67) 

Inserting this is in (|66"D leads again, in the limit kt c ,ujt c — > 0, to the tight coupling 
result. In first order kt c (the square bracket in Eq. ( |6"T| ) has to be expanded up to 
order (kt c ) 5 ) we obtain the dispersion relation 

io = — Z7, with 7 = — ; ~ . (68) 

h ' 20(1 + i?) 3 V ; 

In contrast to the scalar case, vector perturbations show no oscillations (Re(u;) = 0) 
but are just damped. The damping occurs at a slightly larger scale than for scalar 
perturbations, where ■jscaiar — k 2 t c /6 ||PT| . 

The ratio R = 3pt/(4p 7 ) is smaller than ~ 1/4 until the end of recombination. 
We therefore obtain a damping factor / for a given scale k 

(Ik 2 /-WO) \ , 
/~exp — / t c dt) , (69) 



20 Jo 

where t en d(k) is the time at which our approximation kt c < 1 breaks down, i.e., 
kt c (t en d{k)) = 1. The time over which the damping is active is the order of the 
thickness of the last scattering surface, At ~ t dec (Az/ Zd ec ) ~ 0.1t dec . The damping 
scale, the scale at which the exponent in Eq. fl69|) becomes of order unity, is about 

kdamptdec ~ 10. (70) 

The harmonic t corresponding to k damp is £ damp = k damp t ~ 10t /t dec ~ 500. 

After the time t end (k), collisions become unimportant for fluctuations with wave 
number k which then evolve freely, suffering only directional dispersion which induces 
a power law damping oc l/(kAt). Reference || discusses only this second effect. Nu- 
merical experience with scalar perturbations, however, shows that they are typically 
both of similar importance. 
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